Genetics in Medicine Open

Advances in high-throughput sequencing and genetic research have expanded the role of genetics in medicine and society. Population-based screening programs, including neonatal and preconception testing, are increasingly implemented globally, alongside the rise of direct-to-consumer (DTC) genetic testing. The "Genetics and the General Public" Ethics Working Group of the French Federation of Human Genetics (FFGH) assessed knowledge and awareness of genetics within the French population through a nationally representative survey (n=3,013) conducted by the polling firm Ipsos bva. Results indicated that 69% of respondents report an interest in genetics, although their level of knowledge remains limited. Most respondents expressed positive attitudes toward genetics, perceiving it as a major source of hope in healthcare. While a majority indicated willingness to undergo genetic testing for medical purposes, they also reported legitimate concerns regarding the potential results. Despite legal restrictions, 12% reported having ordered a DTC genetic test (5% for genealogical; 5% for medical and 2% for both purposes), and 45% of non-users expressed strong interest in this type of test. Notably, there is a substantial lack of awareness regarding the limitations of these tests and the French legal framework governing their use. These findings highlight critical gaps in public knowledge, emphasizing the need for improved genetic education, including incorporating genetics into school curricula and launching targeted awareness campaigns. These initiatives should help clarify the distinctions between clinically validated genetic tests and DTC genetic testing services, addressing both their benefits and their ethical, legal, and scientific limitations, in order to promote informed decision-making.

IMPORTANCE: Most United States children with neurodevelopmental disorders have not received genetic testing aligned with current guidelines. Integration of genetic counselors into non-genetics departments is a potential strategy to improve uptake, but prevalence and details of integrated care models are unknown. OBJECTIVE: To characterize availability, utilization, and perceived need for genetic counselors across non-genetics departments caring for patients with neurodevelopmental disorders DESIGN: Cross-sectional observational department-level survey SETTING: Child neurology, adult neurology, developmental pediatrics, child psychiatry, and adult psychiatry departments at Intellectual and Developmental Disabilities Research Centers PARTICIPANTS: The survey was distributed to 67 departments across 15 institutions. The departmental response rate was 52% (35/67), with at least one response from 87% (13/15) of institutions. EXPOSURE: Presence/absence of dedicated genetic counselor(s), where "dedicated" was defined as hired by the department MAIN OUTCOME(S) AND MEASURE(S): This was a descriptive study only, with no comparative statistical analyses due to the exploratory nature. RESULTS: One third of departments (34%; 12/35) reported having dedicated clinical genetic counselors. Prevalence was highest in child neurology (67%; 8/12), followed by adult neurology (40%; 2/5) and developmental pediatrics (22%; 2/9), with none in child psychiatry (0/7) or adult psychiatry (0/2). In almost all departments with genetic counselors (92%; 11/12), they directly billed for their services, which universally included pre-test counseling/consent and post-test counseling. In departments without genetic counselors, only 39% (9/23) reported providers ordered their own genetic testing. Among all departments, over half (57%) were interested in adding/increasing genetic counseling support, while 26% were unsure and 17% uninterested. Insufficient funding was the most cited barrier; only one department reported insufficient need. CONCLUSIONS AND RELEVANCE: Though currently implemented in only one third of departments, our findings suggest those with dedicated genetic counselors directly pursue genetic testing (without referring to genetics) more than those without genetic counselors. Interest in increasing or adding genetic counseling support was high, and though funding was a reported barrier, feasible funding models were described. In the context of limited medical geneticists and expanding precision therapies, alternate delivery models for neurodevelopmental genetic testing including genetic counselor integration in non-genetics departments may help to scale and sustain uptake.

Aims Metformin remains the primary treatment for type 2 diabetes, yet over 40% of patients fail to maintain glycaemic control. We aimed to identify patients unlikely to respond to metformin prior to treatment initiation and to evaluate whether on-treatment management can improve glycaemic outcomes in suboptimal responders, informing early treatment decisions. Materials and Methods We analyzed 59,881 longitudinal HbA1c measurements from 7,105 patients with type 2 diabetes receiving metformin monotherapy using real-world electronic health records from Kaiser Permanente Northern California with up to six years of follow-up. We integrated demographic, clinical, genetic, and pharmacological factors to characterize metformin responder phenotypes and quantify the impact of adherence and weight control on time to glycaemic failure. Results Three distinct trajectory-based phenotypes were identified: good (63.6%), poor (8.9%), and non-responders (27.5%). Poor responders initially achieved glycaemic targets but lost control within 2.5 years, while non-responders showed minimal HbA1c reduction and failed within 1 year. Five baseline factors-HbA1c, age at diagnosis, body mass index, sex, and estimated glomerular filtration rate-classified phenotypes with good discrimination (area under the receiver operating characteristic curve = 0.84). Incorporating on-treatment HbA1c further enhanced identification of non-responders. Among suboptimal responders, weight control and improved adherence delayed glycaemic failure by approximately 7 months; however, eventual glycaemic failure remained likely. Conclusions We characterized three clinically relevant metformin responder phenotypes and showed that suboptimal responders can be identified early using baseline features. Poor and non-responders are unlikely to achieve durable glycaemic control with metformin alone and may require alternative treatment strategies.

Purpose: Pathogenic variants in NFIX cause Marshall-Smith syndrome and Malan syndrome (MALNS). We identified a severe subtype of MALNS characterized by adolescent-onset musculoskeletal deterioration and investigated functional consequences of underlying variants. Methods: Clinical data were collected from seven individuals with pathogenic NFIX variants. Wild-type and mutated recombinant NFIX DNA-binding domains (DBDs) were evaluated using biochemical, structural, and DNA-binding assays. Results: Six individuals carrying R116W, R116P, K125E, or G147E NFIX substitutions developed progressive muscle wasting, markedly reduced body mass index, and rapidly progressive scoliosis after the typical childhood features of MALNS; two died from disease-related complications. A seventh individual with R116G did not develop this severe phenotype. Functional studies on recombinant NFIX DBDs showed complete or near-complete loss of DNA-binding activity for R116W, R116P, K125E, and G147E despite preserved protein folding, consistent with disrupted DNA recognition and a potential dominant-negative mechanism. In contrast, R116G exhibited a 7.7{degrees}C decrease in thermal stability, which may support haploinsufficiency mediated by protein degradation. Conclusion: Specific NFIX missense variants define a severe subtype of MALNS associated with progressive musculoskeletal deterioration. In vitro functional studies support variant-specific disruption of DNA binding, providing a mechanistic basis of genotype-phenotype correlations and informing prognosis, clinical surveillance, and therapy development.

Importance Pharmacogenomic (PGx) guidelines can improve medication efficacy and reduce toxicity, but their application in pediatric intensive care units (PICUs) remains largely unexplored. Objective To determine the frequency of medications with established PGx guidelines administered in the PICU and assess the capacity of exome sequencing to capture PGx phenotypes for these medications. Design Retrospective cohort study integrating electronic medical record and exome sequencing data. Setting Morgan Stanley Children's Hospital of NewYork-Presbyterian, a single center tertiary care children's hospital. Participants A total of 4,939 children admitted to the PICU (2020 - 2024), and 192 children admitted to the PICU who underwent exome sequencing for research purposes (2015 - 2023). Exposure Critical illness requiring PICU admission. Main Outcomes and Measures Frequencies of administration of medications with established PGx guidelines in the PICU and the proportion of individuals with exome sequencing with identifiable PGx phenotypes. Results Among 4,939 PICU patients, 37.2% (n=1,837) received at least one medication with established PGx guidelines and 14.4% (n=712) received two or more such medications. Twenty PGx genes were implicated; CYP2C9 was most common (17.3%, n=853). An estimated 8.2% of patients received medications for which PGx-guided recommendations would have altered clinical management. Among 192 patients who underwent exome sequencing, at least one metabolizer phenotype was identified in 62% (n=119). Conclusions and Relevance Many critically ill children receive medications with established PGx guidelines. This study highlights an opportunity for more personalized medicine for critically ill children admitted to a tertiary care hospital and assesses the strengths and weaknesses of exome sequencing to uncover pertinent PGx phenotypes.

Biobank-scale genomic analyses remain computationally expensive, CPU-bound workflows, particularly when adjusting for confounding. Here, we present CuGen, a GPU-accelerated framework for large-scale genomics. CuGen uses UltraLasso, a novel hierarchical application of univariate-guided sparse regression (uniLasso), to select a compact, phenotype-informed active set of fewer than 30,000 variants. This achieves robust leave-one-chromosome-out (LOCO) confounding control, enabling both downstream GWAS and in-sample fine-mapping. Additionally, we introduce the .cugen file format, a genotype representation designed for memory-optimized, high-throughput streaming and random access on GPU hardware. Building on this substrate, we provide a general GPU-accelerated genomics toolkit handling polygenic prediction, data manipulation, quality control, analysis, and visualization. We demonstrate CuGen's efficacy in the UK Biobank with up to 408,624 individuals, where the full GWAS pipeline and fine-mapping against 6.8 million imputed variants completes in approximately 10 minutes on a single high-throughput GPU with 80 GB of memory. The pipeline scales efficiently to massive phenome-wide analyses with sublinear resource consumption.

Summary It has become standard practice to visualize regional signals from genomewide association studies GWAS using LocusZoom plots Similarly GWAS signals are compared to regionally matched quantitative trait loci QTLs ie varianttogene regulation data using LocusCompare plots to aid assessment of candidate traitrelated genes Despite broad usage these tools annotate variants by linkage disequilibrium LD to a single lead or index variant This singleindex representation has limitations for visualizing complex loci that contain multiple independent signals We present LocusBlend an interactive web application for multiindex LDblended visualization of genomic loci LocusBlend supports one or two genomic association summarystatistic datasets and one to three index variants multiindex LocusZoom colorblended plots and matching LocusCompare visualizations Applications to Alzheimers disease GWAS and QTL signals illustrate LocusBlend enables visualization and separation of independent signals despite shared LD and high genomic complexity Overall LocusBlend is aimed at supporting researchers handle the continuously expanding complexity of human genomics findings Availability and Implementation LocusBlend is freely available at httpslocusblendwustledu Publication ready plots are generated in 1min Source code documentation example datasets input templates and reproducibility instructions are available at httpsgithubcomBelloyLabLocusBlend LocusBlend is implemented in Python using Streamlit Plotly and PLINK Supplementary Information Supplementary data are available online

Chronic kidney disease and heart failure disproportionately burden populations of African ancestry, yet Mendelian randomisation (MR) studies of the causal relationship between kidney function and heart failure subtypes have been conducted exclusively in European ancestry populations. We performed a forward two-sample MR analysis to evaluate the causal effect of genetically predicted estimated glomerular filtration rate (eGFR) on heart failure with preserved ejection fraction (HFpEF) and heart failure with reduced ejection fraction (HFrEF) in individuals of African ancestry. Genetic instruments were selected from an African ancestry eGFR genome-wide association study (N = 67,943) at genome-wide significance, with linkage disequilibrium clumping using an African ancestry reference panel. Heart failure subtype summary statistics were obtained from the Million Veteran Program (HFpEF: 5,379 cases / 113,041 controls; HFrEF: 9,104 cases / 109,632 controls). Six independent SNPs (F-statistics 30.5 &#8211 107.3; R&#178 = 0.62%) were retained as instruments. The primary inverse-variance weighted analysis provided no evidence of a causal effect of eGFR on HFpEF (OR 0.92, 95% CI 0.80 &#8211 1.06, p = 0.248) or HFrEF (OR 0.98, 95% CI 0.78 &#8211 1.23, p = 0.878). Sensitivity analyses were directionally consistent. There was no evidence of heterogeneity or directional pleiotropy. Minimum detectable effects at 80% power were OR 1.28 for HFpEF and OR 1.22 for HFrEF. These null findings should be interpreted as inconclusive given current power constraints; larger ancestry-matched studies are needed.

Background: Actigraphy-derived rest-activity rhythm (RAR) features are widely used to characterize clinical states in bipolar disorder (BD). Both mean levels and temporal variability of these features have been associated with mood episodes; however, variability measures are often statistically coupled with the mean, particularly in skewed distributions. This raises a question as to whether variability reflects a separate characteristic of the data or whether the observed association arises from statistical properties of the data. Objective: In this study, we aim to determine whether temporal variability of actigraphy-derived RAR features provides standalone information on mood episodes in BD beyond mean activity levels after accounting for mean-variance dependence. Methods: We analyzed actigraphy data from a subset of 72 participants with BD drawn from a larger longitudinal study, extracting 22 daily RAR features aggregated weekly as sample mean (MEAN) and within-week temporal variability computed as sample standard deviation (VAR). Variance-stabilizing transformations (Box-Cox or Yeo-Johnson) were applied to the entire study cohort to reduce mean-variance dependence. Associations with mood episodes and remission (mania: n=34; depression: n=58 annotated participants) were evaluated using generalized linear mixed-effects models with a logistic link function, including univariate (MEAN or VAR) and multivariate (MEAN+VAR) specifications, assessed by likelihood-based metrics and the area under the receiver operating characteristic curve (AUC). Results: Transformations reduced mean-absolute correlations from 0.43 to below 0.06. Temporal variability remained significantly associated with clinical state for 11/22 RAR features in mania and 16/22 features in depression, with all significant associations remaining after false discovery rate correction (p<0.05). Joint models showed modest incremental gains (AUC 3%-4% overall; up to 12% in mania, 7% in depression), with absolute performance remaining limited (AUC 0.50-0.66). In both mania and depression, nearly all significant variability-based regressors contributed incremental information beyond mean-based models. Only sleep duration and activity changes around wake time (+-1 hour), did not improve discrimination between mania and remission. Conclusions: Temporal variability in RAR features can be considered a standalone state marker of mood episodes not captured by mean activity. We found it to be more consistently associated with depression than mania. Its incremental discriminative contribution is modest, suggesting greater utility within multivariate or multimodal frameworks.

We present FootNet, a 453-image multi-view smartphone foot dataset for binary foot segmentation, with expertannotated masks across six anatomical views (dorsal, medial, and plantar, both left and right). We benchmark four segmentation models under a controlled protocol: U-Net with a MobileNetV2 encoder achieves the best performance (IoU 0.9268, Dice 0.9608, 95 % CI [0.9209, 0.9320]); DeepLabV3 with MobileNetV3-Large scores IoU 0.8984 (Dice 0.9449); UNet++ with MobileNetV2 scores IoU 0.8913 (Dice 0.9391); and SAM ViT-B with oracle boundingbox prompt scores IoU 0.9219 on the matched 191-image subset. Bonferroni-corrected Wilcoxon signed-rank tests (k = 6 comparisons) show U-Net significantly outperforms DeepLab (p < 0.001, r = 0.638) and SAM ViT-B with oracle boundingbox (p = 0.005, r = 0.202); UNet++ does not significantly differ from DeepLab (p = 0.062). Connected-component postprocessing yields negligible benefit (mean {triangleup}IoU = +0.0003, 12 of 453 images improved). The extended dataset is available upon request

Lithium has been the gold standard for the treatment and prevention of relapse in bipolar disorder for over 60 years. Guidance from the National Institute for Health and Clinical Excellence states explicitly to 'offer lithium as a first-line, long-term pharmacological treatment for bipolar disorder'. Yet, in the last two decades its use has been in decline with clinicians favouring anticonvulsants or antipsychotics when treating this condition. In this study, we have used three openly available datasets containing prescribing data from primary and secondary care to explore trends in the use of lithium in England, showing both regional and temporal variance between 2015-2024. We have shown that lithium use declined in primary care by 20.9% in the last ten years (2015-2024) and 10.9% overall in the last five years (2019 to 2025). We have also shown how there is some regional variation in the source of lithium for patients, although the vast majority is prescribed in primary care. Further research into clinical behaviour is needed to understand what is driving the decrease in lithium usage, and what barriers and enablers may influence its use across the country.

Background. When a GLP-1 patient stops responding, should the clinician push the dose or change the drug? Observational answers conflate two distinct sources of confounding. Most early-week "escalations" in real-world data are FDA-mandated titration steps rather than deliberate clinical decisions, and patients who deviate do so for reasons we cannot observe. Semaglutide and tirzepatide are also not equivalent milligram-for-milligram, so naive class-switch comparisons mix mechanism and dose. We resolve both by restricting to post-titration patients and reclassifying treatments under the Whitley 2023 dose-equivalence framework. Methods. From 68,969 telehealth GLP-1 patients we built a post-titration cohort. Each patient's index time is the day they completed at least four weeks at therapeutic dose (Whitley tier 3 or higher: semaglutide 1.0 mg or tirzepatide 5 mg). Confirmed slow response is less than 5% total weight loss at the index, consistent with FDA weight-management drug-development guidance and AACE/ACE criteria. We compared four post-index strategies against continuing the current regimen: within-class dose escalation, equipotent class switch (a Whitley tier change of 1 or fewer), and class switch with potency increase. Direction-specific analyses split switches into semaglutide-to-tirzepatide and tirzepatide-to-semaglutide arms. Outcomes were percent weight loss at 12 and 24 weeks post-index. We estimated effects six ways: propensity-score matching; IPTW with linear and gradient-boosted propensities; the g-formula with linear and gradient-boosted outcome models; and AIPW, the doubly-robust estimator we use as the tiebreaker. Two-layer inverse probability of censoring weighting addressed strategy adherence and outcome ascertainment. We computed E-values, ran a negative-control specification, and stratified by tolerability. Results. The post-titration cohort comprised 24,876 confirmed slow responders. Within-class dose escalation produced a small consistent benefit at 24 weeks: AIPW +0.64 pp (95% CI +0.16 to +1.12), with five non-AIPW estimators ranging +0.47 to +0.76 pp. The continue arm itself lost an additional 8.07 pp over the same window (96% continued to lose), so escalation is a marginal addition to a substantial natural slope, not a rescue. Equipotent class switching from semaglutide to tirzepatide was inconclusive: linear and matching estimators ranged +1.26 to +1.80 pp, but AIPW was -0.33 pp (95% CI -1.27 to +0.60) with only 90 treated patients and limited propensity-score overlap. Class switch with simultaneous potency increase (sema to tirz) gave AIPW +0.65 pp at 12 weeks (95% CI +0.37 to +0.92, n = 80). A negative-control specification yielded ATE -0.13 pp, indicating the pipeline did not generate spurious signal. A held-out-fold prognostic-threshold sensitivity gave a null effect (+0.04 pp), correcting an earlier circular +1.06 pp estimate. Conclusions. Among confirmed post-titration slow responders, within-class dose escalation adds approximately 0.6 percentage points at 24 weeks on top of an 8 percentage point natural slope, consistently across six estimators including doubly-robust inference. This headline effect is small and not robust to modest unmeasured confounding (E-value 1.27) or to MNAR-style outcome attrition (tipping point delta approximately 1.2 pp), so it should be read as hypothesis-generating rather than practice-changing. Class-switching evidence is inconclusive; linear-estimator results suggesting benefit did not survive doubly-robust estimation in small treated samples with limited propensity overlap. The dose-ladder framework, with phase-specific evidence grading, is hypothesis-generating and insufficient on its own to change practice.

Objective: There is a bidirectional relationship between emotional distress and pain. However, this relationship is understudied in people with HIV in low-resource settings. This study sought to describe the temporal relationship between emotional distress and pain in people with HIV. Design: Longitudinal observational study. Methods: Participants with virally suppressed HIV, reporting either no pain or persistent pain at baseline, provided weekly remote ratings of distress, worst pain, and average pain using 0-10 visual analogue scales. Within-individual fluctuations in distress and pain were visualised over time. Group-level correlations were determined using Spearman's correlation tests. Cumulative link mixed models assessed whether distress and pain each predicted the other in the following week. Results: 72 participants provided responses over 49 weeks. The participants had a median (IQR) age of 43 (37-51) years, 63% (n=45) were unemployed and most were females (n=51;71%). Distress and pain fluctuated concurrently within individuals: distress was positively correlated with worst pain ({rho}=0.66, 95% CI= 0.60-0.72, p<0.001) and average pain ({rho}=0.70, 95% CI=0.64-0.75, p<0.001) intensity within the same week. Worst pain (OR=1.42, 95% CI=1.17-1.71, p<0.001) and average pain (OR=1.43, 95% CI=1.20-1.71, p<0.001) intensity both predicted distress in the next week. Distress predicted worst pain intensity (OR=1.25, 95% CI=1.07-1.46, p=0.023) but not average pain intensity (OR=1.19, 95% CI=1.01-1.40, p=0.152) in the next week. Conclusions: The temporal relationship between distress and worst pain intensity was bidirectional, whereas distress did not temporally predict average pain intensity. Both pain and emotional distress should receive attention from HIV research and clinical care in low-resource settings.

Background Post-infection conditions (PICs), such as Long Covid, are associated with heterogeneous, fluctuating symptoms that profoundly affect daily functioning. Despite moderate-certainty evidence from the NIHR-funded LISTEN trial (COV-LT2-0009) that personalised self management support improves outcomes and may reduce societal and economic impacts of Long Covid, many people living with PICs still receive condition-specific services, generic advice, or stand-alone digital tools that do not address their complex needs. Aim To map care approaches in general practice and synthesise UK evidence for PIC management. Design and setting Scoping review and online survey. Method A two-phase study was conducted: (1) a scoping review of UK evidence on PIC management in general practice; and (2) a supplementary online survey of practitioners working in UK general practice to provide contextual insights. Results The scoping review identified 32 studies focused on Long Covid. One study included a comparator group (ME/CFS). Study populations were predominantly white ethnicity and female. Evidence for non-Covid PICs in UK general practice was largely absent. The supplementary survey (n=46) provided preliminary practice-level insights. Healthcare practitioners reported varied PIC presentations, diagnostic uncertainty, limited referral pathways, inequitable access, and low confidence in managing PICs. Conclusion Evidence informing PIC management in UK general practice remains predominantly Long Covid-focused and may not reflect the range of PICs encountered in practice. While survey findings are preliminary and require confirmation in larger samples, they highlight uncertainty around PIC management. Further research is needed to evaluate whether existing Long Covid pathways should be expanded or complemented by broader PIC models. Keywords general practice; Long Covid; self-management; post-viral syndromes

Background: Children born prematurely (before 37 weeks) or admitted to the neonatal unit (NNU) are at increased risk of adverse long-term physical health outcomes. It is also recognised that there is an association with later academic performance and special educational needs, however it is not clear whether these broad risk factors could be used as stand-alone heuristics to identify children who may benefit from additional support in educational settings. We aimed to examine the associations between neonatal unit (NNU) admission and educational attainment in mid-childhood. Methods and Findings: Pregnancy data from a prospective birth cohort (Pregnancy Outcome Prediction Study, Cambridge, United Kingdom, 2008-2012) were linked to national educational outcomes (Department for Education, United Kingdom). Multivariable regression models adjusted for maternal, child, and socioeconomic factors were used to evaluate associations between (i) all NNU admissions, (ii) at term NNU admissions >48 hours, (iii) preterm birth without ongoing physical health needs, and educational outcomes at ages 5-11 years. Children who required any NNU care were more likely not to meet expected educational standards across multiple ages and domains in early and mid-childhood: age 5 early year foundation (aOR 1.64, 95% CI 1.19-2.27, p=0.003), phonics at age 6 (aOR 2.43, 95% CI 1.72-3.57, p<0.001), and at age 7 (here assessments were divided into multiple domains): reading (aOR 1.67, 95% CI 1.18-2.38, p=0.004), writing (aOR 1.72, 95% CI 1.25-2.38, p<0.001), mathematics (aOR 1.56, 95% CI 1.09-2.22, p=0.020), and science (aOR 1.85, 95% CI 1.22-2.78, p=0.003). Similar patterns were observed among both at term-born infants who stayed >48hrs in NNU (phonics assessment at age 6 aOR 2.26, 95% CI 1.51-3.36, p<0.001) and in children born preterm without long-term physical health sequelae (phonics assessment at age 6 aOR 3.07, 95% CI 1.96-4.81, p<0.001). These associations were robust to adjustment for demographic, perinatal, and socio-economic factors. By age 11, differences in academic attainment were attenuated and no longer clearly distinguishable across all exposure groups. However, there was an increased likelihood of special educational needs (SEN) at age 11 associated with any NNU admission (aOR 1.78, 95% CI 1.15-2.73, p=0.009), at term NNU admission for >48hrs (aOR 1.88, 95% CI 1.19-3.00, p=0.007), and children born preterm without long-term physical health sequelae (aOR 1.50, 95% CI 1.00-2.25, p=0.049). Predictive performance of any NNU admission for SEN at age 11 was moderate (AUC 0.70, 95% CI: 1.14-2.65, p=0.010), with balanced sensitivity and specificity and high negative predictive value. Conclusions: NNU admission, for both term and preterm infants, is associated with poorer educational outcomes and an increased likelihood of special educational needs in mid-childhood.

Estimating dengue force of infection (FOI) is essential for understanding transmission dynamics and targeting intervention programmes, yet surveillance data in endemic settings required for estimations are often incomplete, with varying formats. We developed a Bayesian hierarchical catalytic model that jointly fits age-stratified case data, aggregate case data, and seroprevalence surveys within a single framework, incorporating external covariates to improve parameter identifiability. Synthetic validation showed that covariates alone recovered accurate FOI point estimates even when most districts contributed only aggregate data, but did so with poorly calibrated uncertainty; anchoring the model with a single seroprevalence survey was necessary to bring credible interval coverage close to nominal. Applied to 128 districts across Java and Bali, Indonesia (2016-2024), the model revealed substantial spatial heterogeneity in FOI and reporting rates. Many districts in Java exceeded the WHO-suggested seroprevalence threshold for vaccine introduction, yet were classified as low-priority when using reported incidence as prioritisation criterion, particularly in areas with weak surveillance. Model-based seroprevalence estimation, integrating multiple data sources, offers a more consistent basis for identifying high-priority districts for vaccine introduction, and is less susceptible to surveillance bias than reported incidence.

Stochastic epidemic models are a cornerstone of infectious disease epidemiology and are often used to study intervention scenarios. However, large run-to-run variability can make intervention effects difficult to estimate precisely. We revisit the epidemic Sellke construction, which assigns each individual an infection threshold for the cumulative infection hazard such that, conditional on the thresholds, the epidemic trajectory becomes deterministic. This enables coupling of simulations with and without an intervention, yielding low-variance effect estimates even when outcomes such as final size or peak incidence vary widely between runs. We develop an exact, event-driven implementation that maintains infection and recovery events in priority queues. Cumulative infection-hazard updates require O(log N) time per event, yielding overall complexity O(Elog N) for E events in a population of size N. The implementation achieves computational performance comparable to the classical Gillespie algorithm while naturally accommodating non-Markovian infectious periods and complex infectiousness profiles. We illustrate the approach using distance-dependent spread of avian influenza between poultry farms in the Netherlands and a multilayer population with households, schools, and workplaces. In both examples, coupling enables efficient within-run comparisons of intervention scenarios across stochastic realisations.

Background: Opioid use disorder (OUD) is heritable, yet most genome-wide association studies (GWAS) have focused on European populations, leaving the genetic architecture of OUD in non-European populations underexplored. Methods: We conducted GWAS of OUD across three ancestries using electronic health records and genomic data from 52,357 All of Us Research Program participants (8,912 cases; 43,445 matched opioid-exposed controls; 48.5% female). Participants were stratified into European (EUR), African (AFR), and Admixed American (AMR) ancestry groups for logistic regression GWAS, with independent replication in the Million Veteran Program. We then applied the deep-learning model AlphaGenome to predict the tissue-specific transcriptomic and splicing consequences of top risk variants across 13 reward-pathway brain regions. Results: We identified and replicated a novel DDX6 risk locus, alongside established OPRM1 and FURIN signals. AlphaGenome predicted the DDX6 regulatory allele downregulates the stress-resistance gene FOXR1 in the nucleus accumbens, while the protective OPRM1 variant (rs1799971) upregulates OPRM1 expression across reward networks. Other signals of interest included IL6R and SHISA9 (EUR); GHR (AFR); and ASTN2 (AMR). Conclusions: This study identifies DDX6 as a novel OUD risk locus, replicates associations with OPRM1 and FURIN, and highlights biologically plausible ancestry-specific signals in AFR and AMR populations. We also replicated top variants in an independent population. Finally, integrating GWAS with deep-learning annotations provides specific, localized biological hypotheses to guide future experimental validation and targeted therapeutics.

Background: Monoclonal antibodies have emerged as a prophylactic strategy to prevent symptomatic SARS-CoV-2 infection in immunocompromised individuals. However, the evolutionary and clinical implications of breakthrough infections under this regime remain unclear. Methods: A male in their 80s with a haematological/oncological diagnosis received a 2000 mg intravenous infusion of sotrovimab in March 2023 and was diagnosed with COVID-19 by RT-qPCR from a nasopharyngeal swab in August 2023. Weekly samples (n=24) were collected through February 2024 (171 days). All samples underwent whole-genome sequencing, with select mutations subjected to functional assessment. Findings: Sequencing identified the GE.1 lineage at all timepoints. An intra-host recombination event in ORF1ab (positions 8942-12458) was detected prior to 23 weeks post-detection, followed by a 14-fold increase in viral load (7.42e+06 to 1.00e+08 RNA copies/mL) and a marked shift in the viral population. E340D, a sotrovimab resistance mutation, was detected at low abundance (46%) within the first week post-infection, fluctuated over time, and was nearly fixed by week 15 (107 days) post-detection. We assessed five spike mutations - V36M, S98F, and V213G in the N-terminal domain, Y505P in the receptor-binding domain, and P681Q near the S1/S2 cleavage site - and additionally evaluated the impact of E340D. V36M conferred the highest infectivity across all cell lines, with the most significant effect in low-TMPRSS2 cells. While all mutations showed enhanced infectivity with the addition of E340D, the effect was most pronounced in mutations with lower baseline infectivity. The addition of E340D significantly decreased relative neutralizing titres for V36M, S98F, and V213G, enabling escape from neutralizing antibodies in XBB-responsive individuals, illustrating an enhanced phenotypic advantage. Patient neutralizing activity was absent pre-sotrovimab, and sotrovimab-induced neutralization was further compromised by selection of E340D. Interpretation: Sotrovimab pre-exposure prophylaxis in an immunocompromised patient did not prevent SARS-CoV-2 infection, and selected for resistant mutation E340D, with unexpected fitness consequences across non-receptor binding domain spike regions.

Insulin signalling contributes to neurodevelopment and brain function, and insulin resistance (IR)-related traits are associated with cognitive performance. However, the genetic architecture shared across specific cognitive domains and IR-related phenotypes remains insufficiently defined. We analysed large-scale genome-wide association study summary statistics for 11 IR-related traits (N=53,334-933,970) and 10 cognitive measures (N=28,156-436,853) to quantify global and local genetic correlations, fine-map shared association signals, and annotate implicated genes and drug-gene interactions. Pairwise global and local genetic correlations were estimated, and shared high-confidence variants were prioritised using the multivariate Sum of Single Effects model. Positional and expression quantitative trait locus mapping was performed, and implicated genes were examined through functional annotation, tissue enrichment, and drug-gene interaction analyses. Low-to-moderate genetic correlations were observed between six IR-related traits and seven cognitive measures (|rg|=0.08-0.34), with predominantly opposite directions, except for correlations involving visual declarative short-term memory. Local genetic correlations showed mixed effect directions across most trait pairs, and multivariate fine-mapping prioritised 696 shared likely causal variants with high posterior support. Gene annotation indicated enrichment in several pathways, including immune-related, signal transduction, neurogenesis, neurotransmitter metabolism, receptor regulation, and lipid and cholesterol metabolism regulation. Implicated genes were expressed across various brain regions and showed prior associations with neuropsychiatric and cardiometabolic conditions. Several drug-gene interactions were identified, involving immunomodulatory and anti-inflammatory compounds. These findings indicate widespread heterogeneous genetic overlap between IR-related traits, particularly body mass index and waist-to-hip ratio, and cognitive measures of general intelligence, processing speed, and short-term visual declarative memory. The findings prioritise apolipoprotein-related lipid transport and inflammatory and oxidative stress pathways as candidate mechanisms linking cognitive, cardiometabolic, and neuropsychiatric phenotypes.